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ABSTRACT 

By using a decomposition of the transfer matrix of the (/-state Potts Model on a 
three dimensional m x n x n simple cubic lattice its determinant is calculated exactly. By 
using the calculated determinants a formula is conjectured which approximates the critical 
temperature for a d-dimensional hypercubic lattice. 
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I. Introduction 



The two dimensional (/-state Potts models [1,2] for various q have been of interest as 
examples of different universality classes for phase transitions and, for q = 3,4 as models 
for the adsorption of gases on certain substrates [3,4,5]. For q > 3 the free energy has never 
been calculated in closed form for arbitrary temperature. Some exact results have been 
established for the model: from a duality relation, the critical point has been identified 
[1]. The free energy and latent heat [6,7,8], and magnetization [9] have been calculated 
exactly by Baxter at this critical point, establishing that the model has a continuous, 
second order transition for q < 4 and a first order transition for q > 5. Baxter has also 
shown that although the q = 3 model has no phase with antiferromagnetic long-range order 
at any finite temperature there is an antiferromagnetic critical point at T = [9]. The 
values of the critical exponents (for the range of q where the transition is continuous) have 
been determined [10-12]. Further insight into the critical behaviour was gained using the 
methods of conformal field theory [13]. Our knowledge about the three dimensional Potts 
model is much less than the two dimensional case. The three dimensional 3-state Potts 
ferromagnet serves as an important model in both condensed matter as well as high energy 
physics [2]. Experimental realizations are structural in some crystals, and theoretically this 
model has attracted much interest as a simple effective model of finite-temperature pure- 
gauge QCD. Consequently it has been studied in the past few years by many authors using 
quite a variety of different techniques, for an incomplete list see [14-20, and refs therein]. 
By using Monte Carlo simulations with the help of standard finite size scaling methods 
the characteristic parameters of the phase transition (transition temperature, latent heat, 
etc) have been estimated with varying accuracy [21]. As a result there is by now general 
consensus that this model undergoes on a simple cubic lattice a weak first order phase 
transition from a three fold degenerate ordered low temperature phase to a disordered 
phase at high temperature. However, our knowledge about exact value of the partition 
function and spectrum of the three dimensional q-state Potts model is almost nothing. It 
is thus of continuing value to obtain further information about the three dimensional Potts 
model. In this paper by using edge transfer matrices the product of the eigenvalues for the 
transfer matrices of the (/-state Potts model on the three dimensional lattices is calculated. 
The paper is organized as follows: In section II, by using the standard representation 
of the transfer matrices for the (/-state Potts model on 2 x n x n lattices and using a 
method which is suggested by R.J.Baxter [29], the determinants of the transfer matrices 
are calculated exactly and then the calculation is generalized to the three dimensional 
m x n x n simple cubic lattices. Finally from the general relation for the determinant a 
formula for the <i-dimensional (/-state Potts model is conjectured which approximates the 
critical temperature with a few percent error. 



II. Determinant of the Potts model transfer matrix 



The (/-state Potts model has served as a valuable model in the study of phase transition 
and critical phenomena. On a lattice, or more generally on a graph G, at temperature T 
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this model is defined by the partition function: 

Z(G,q 1 k)=J2 e ~' 3H (!) 

with the Hamiltonian 

H = -J Sfaaj) (2) 

<i,j> 

where 8(si,sj) is the kronecker delta and Sj = 1, . . . , q are the spin variables on each vertex 
% G G , (3 = (^T)" 1 , k = (3J; and < i,j > denotes pairs of adjacent vertices. Consider 
an n x n square lattice with periodic boundary conditions. For the two dimensional Ising 
model, partition function can be written as trace of the product of transfer matrices and 
the eigenvalues can be calculated exactly [22-27]. There are several representations for the 
transfer matrix of the two dimensional g-state Potts model [22,28] which can be used to 
represent the transfer matrices for m-layer models. Determinant of the transfer matrix for 
the two dimensional Potts model can be calculated by different methods [29,30,31]. It is 
quite hard if not impossible to generalize our method [30] to higher dimensions (d > 2 
). We will use a method which is suggested by R.J.Baxter and generalize it to calculate 
the determinant of the transfer matrix for the q-state Potts model on three dimensional 
m x n x n simple cubic lattices. At first we calculate the determinant for the two-layer 
lattices, generalization to m-layers is simple and straightforward. For a two layer lattice the 
transfer matrix element connecting the stripes s±, . . . , s n , t±, . . . , t n and s[, . . . , s' n , t[, . . . , t' n 
can be written as 

< Si, . . . , s n , £i, . . . , t n | P | s^, . . . , s n , . . . , t n >= 

n 

J~J e k[8(s k ,s k + 1 )+S(t k ,t k + 1 )+S(s k ,t k )] e k[S(s k ,s' k )+S(t k ,t' k )] ^ 
k=l 

We consider the stripe s±, . . . , s n , t±, . . . , t n in a plane with a periodic boundary condition 
along one direction (x) and free boundary conditions in the other direction (?/) , i.e. s n j r \ — 
si,t n+ i = t\. The transfer matrix operates along the z direction which is perpendicular 
to the x — y plane and connects the two stripes. Let us define two q 2n x q 2n matrices V2 
and V{ whose matrix elements are given by 



<s 1 ,...,t n \V 2 \s' 1 ,...,t' n >= S a , a > S ttt > J] e ms k ,s k+1 )+8(t k ,t k+1 )+8(s k ,t k )] (4) 

k=i 

n 

< Sl ,...,t n I V{ I s[,...,t' n >= Y[ e k ^ s ^+ s ^ t '^ (5) 

k=i 

where 
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1=1 1=1 

Both V 2 and V/ are products of n matrices [29] : 

v 2 = w 12 w 23 ...w nA 
v{ = u 1 u 2 ...u n 

where W v ^+i is a diagonal matrix with entries 

[W- M+1 ] sM ^ = e M«(«i,« i+ i)+«(t i ,ti+i)+«(- i ,ti)] < j flja , ( y tit/ 
and Ui (i = 1, . . . , n) is the matrix with entries 

n 



(6) 



(7) 



(8) 



(9) 



(10) 



Simultaneously permuting si, . . . , £ n and s[, . . . , t' n is equivalent to a rearrangement of the 
rows and columns of these matrices, so 



det W 12 = det W 23 = ■ ■ ■ = det W ni 



(11) 



det Ui = det Z7 2 = • • • = det l/ n 



(12) 



Hence 



det P = (det W 12 det ?7i) n (13) 
Wi 2 is a diagonal matrix and its determinant is the product of its diagonal entries 



det^ 12 = n n 



3 /e5(5i,52)«fe5(ti,t2)«fe5(si,*l) 



s lv) S n ^ 1 ? ■ ■ 



nn 

Si,S 2 tl,*2 

c 2n-4 

,3fc<? 3 



,fc<5(si,s 2 ) „fe(5(ti,t 2 ) /3 fe(5(si,ti) 



= e 3fc * 



C/i is a block-diagonal matrix, consisting of g diagonal blocks, each of the form 

A = (e Iqxq ~\~ Pqxq) <S> (e Iqxq ~\~ &qXq) 
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(14) 

(15) 

(16) 
(17) 

(18) 



where a is a q x q matrix with zero diagonal elements and unit elements on all other entries 
and / is a q x q unit matrix. The determinant of the q 2 x q 2 matrix A is given by, 



det A = (e k + q - l) 2q (e fc - l) 2( ^-i) 
by using the following definitions 



x 



e k + q — 1 



, y = e 



e k - 1 
we arrive at 

dettfi = (detA) q2n ~ 2 = [x 2q (y-l) 2q2 ] q2 
so the determinant of the transfer matrix is given by 

detP=[x(y-l) q ] 2nq2n ~ 1 y Snq2n ~ 1 



(19) 



(20) 



(21) 



(22) 



For m-layer lattices after a straightforward calculation we get to the following results. The 
matrix A is a direct product of m factors 



A — (e Iqxq + PqXq) <E> (e Iqxq ~\~ O'qXq) 

and so we arrive at 

det A = x mqm ~\y - 1) 



(c -^gXg + CqXq) 



mq 



(23) 



(24) 



m(n— 1) 



(dett/O 71 = (detA) n,? (25) 
is again a diagonal matrix and its determinant is the product of its diagonal entries 



detW 12 = \J 



„(1) „(1) ,("») 



s l >*2 1 ' 2 



m(n-2) 



= (2m-l)fcqr 2 



= e (2m-l)feq r 



— 2) 



(26) 

(27) 

(28) 
(29) 



Finally we arrive at the following formula for the determinant of the transfer matrix for 
the g-state Potts model on a three dimensional m x n x n simple cubic lattice 



[mn — 1] 



detP = [x(y - l)Q]mnq—' y n(2 m -l) q 



[mn — 1] 



(30) 
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Considering a periodic boundary condition in both x and y directions (a plane which is 
perpendicular to the direction that transfer matrix operates) , a straightforward calculation 
shows that the determinant is given by the following simple equation 

det P = [xy 2 (y - l)?]™ 2 ^ 2 - 1 (31) 

There is a temperature where the product of the eigenvalues does not depend on the size 
of lattices. It should be noted that all of the eigenvalues are positive and det P = 1 can be 
simplified to the following equation 

xy 2 (y- 1)^ = 1 (32) 

It is clear that (29) has a solution for k = (3 J which is a function of q and does not depend 
on m or the size of lattices. There is another interesting point about the periodic boundary 
conditions in all directions. Consider a one dimensional lattice with a periodic boundary 
condition, the transfer matrix is given by 

P^ ^ = e Iqxq + CqXq (33) 

The determinant can be calculated easily 

detP (1) =x{y-l) q (34) 

For a two dimensional m x m lattice with a periodic boundary condition the determinant 
is given by [29,30,32] 

detP {2) = [xy{y-l) q ] mqm ~ 1 (35) 

and the determinant for a three dimensional m x m x m lattice with a periodic boundary 
condition for the two dimensional plane which is perpendicular to the direction that the 
transfer matrix operates is given by 

det = [xy 2 (y - l)?]- 2 ^ 2 " 1 (36) 

This calculation can be extended to the d dimensional lattices. Considering a toroidal 
boundary condition for the d — 1 dimensional hypersurface which is perpendicular to the 
direction that the transfer matrix operates we arrive at the following relations 

(detW 12 ) m = y 8 ^- 1 ^' 1 (37) 
(det[/ 1 )- = (det^)^ (s/m)(m " 1) (38) 

det ,4 = x (-/™tf" m) -\ y - (39) 

where s = m d ~ 1 . Finally we arrive at the following general formula for the determinant 
at any dimension (d = 1,2,3,...) for x • • • x hypercubic lattices 
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det P (d) = [xy d -\y - l) q ] sqB ' (40) 

For a two dimensional lattice the critical point corresponds to the maximum of the following 
factor of the determinant. 

-j^w = xy (41) 

We may expect that the maximum of the following factor of the d-dimensional determinant 
will give us an approximate formula for the critical temperature (we cannot prove this 
except for the two dimensional lattices). 

xy d ~ 1 (y — l) q rf , 

—. I = X yd-1 42 

The maximum occurs at the following point 

. , . , r 2(q-l)-d(q-2) + V4q(d-l) + q2(d-2)2 

k m ax{q, d) = In [ 2(d- 1) ^ ^ ' 

So we may conjecture that (43) is an approximation for the critical temperature. Interested 
reader may compare (43) with data which is given in [2] (see tables II and III in page 256) 
and see that it is indeed a good approximation for the critical temperatures. 



(q,d) (44) 

It may be interesting to extend these results to other lattices with different boundary 
conditions. This decomposition of the transfer matrix may be useful for obtaining other 
exact results for the m-layer or three dimensional lattices. 



V. Conclusion 

In this work determinant of the transfer matrix for the g-state Potts model on a three 
dimensional mxnxn simple cubic lattice is calculated exactly and the critical temperature 
is approximated by a conjectured formula. 
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